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Theoretical aspects of potential astrophysical sources of the highest energy 
cosmic rays are discussed, including their energy budget and some issues on 
particle escape and propagation. After briefly addressing AGN jets and GRBs, 
we highlight the possibility of heavy nuclei originating from cluster accretion 
shocks. The importance of X-ray and gamma-ray signatures in addition to 
neutrinos as diagnostic tools for source identification is emphasized. 



1. Introduction 

Several decades after their discovery, the origin of ultra-high energy cosmic 
rays (UHECRs), cosmic particles with energies 10 18 -10 20 eV and above, 
remains one of the biggest mysteries in physics and astrophysics. 1 ' 2 Many 
issues contribute to the difficulty of the problem. For experimentalists, the 
extremely low event rates necessitate detector facilities with huge effective 
area in order to obtain reliable results. The relevant energies far exceed 
those of terrestrial experiments, often making the determination of basic 
observables such as particle energy and composition dependent on inter- 
action models with large uncertainties. On the theoretical side, conceiving 
viable explanations for the production of UHECRs with conventional phys- 
ical mechanisms in known astrophysical objects is a great challenge. The 
unavoidable yet uncertain influence of Galactic and extragalactic magnetic 
fields on UHECR propagation pose further complications. 

However, great advances are expected in the coming years with the ad- 
vent of new generation facilities such as the Pierre Auger Observatory and 
the Telescope Array, as well as future projects such as the Extreme Uni- 
verse Space Observatory. Combined with crucial complementary informa- 
tion from neutrino, X-ray and gamma-ray observatories, the solution of the 
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mystery could be within sight soon. This article, by no means a thorough 
review, discusses selected theoretical topics in this exciting field, focusing 
on the astrophysical aspects. 

2. Issues on UHECR propagation 

We first touch upon some issues concerning the propagation of UHECRs, 
the basics of which have been well reviewed elsewhere (e.g. Ref. 3). The 
observed global isotropy in the arrival directions strongly suggests that 
UHECRs are of extragalactic origin. If UHECRs are protons, photopion in- 
teractions with cosmic microwave background (CMB) photons must induce 
severe energy losses at > 7 x 10 19 cV for propagation distances > 30 Mpc. 
Unless the sources lie much nearer, a spectral ("GZK" 4 ) cutoff is expected 
above these energies. Whether or not there is actual evidence for this in the 
observed UHECR spectrum is a matter of controversy at the moment, and 
will not be discussed here. 

Here, wc call attention to the possibility that the highest energy UHE- 
CRs are composed mainly of heavy nuclei such as iron. For nuclei, the 
dominant energy loss process above 10 19 eV during intergalactic propaga- 
tion is photodisintegration and pair production interactions with photons 
of the far-infrared background (FIRB) and the CMB. 5 Evaluations based 
on recent determinations of the FIRB show that the energy loss distance 
for iron nuclei at 10 20 cV is > 100 Mpc, somewhat larger than that for 
protons. 6 Observationally, the nuclear composition of UHECRs, particu- 
larly at the highest energies, is quite uncertain. Important clues come from 
fluorescence measurements of shower elongation rates by the HiRes experi- 
ment, indicating that the composition changes from heavy-dominated below 
10 18 eV (presumably of Galactic origin) to light-dominated above this en- 
ergy. 7 However, above ~3x 10 19 eV the statistics run out and there is no 
reliable information currently available. Even for the lower energy range, 
other methods give rather different results, 8 and systematic uncertainties 
in interaction models and atmospheric optical attenuation remain a seri- 
ous concern. 9 In any case, at present, it is quite viable that UHECRs at 
> 3 x 10 19 eV are predominantly heavy nuclei originating from extragalactic 
sources. This picture will be elaborated on in Sec. 4. 

Another crucial issue that is often underestimated is deflection by ex- 
tragalactic and Galactic magnetic fields, which can affect the arrival di- 
rections of UHECRs and also act to lengthen their effective propagation 
distance. The strength and distribution of magnetic fields in the intergalac- 
tic medium and at high latitudes in the Galaxy are very poorly known, both 
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observationally and theoretically. Faraday rotation measurements of distant 
radio sources give only upper limits in the nanogauss range for intergalactic 
fields on average, subject to assumptions on the field coherence scale and 
ionized gas distribution. 10 Realistically, whatever their origin, intergalac- 
tic fields can be expected to have some correlation with the distribution 
of large scale structure. Attempts to model this through cosmological sim- 
ulations haven given different results depending on the input physics and 
numerical methods. 11,12 The effect of Galactic magnetic fields is also model- 
dependent. 13 Definitive answers may not come until the operation of the 
Square Kilometer Array (SKA), slated to undertake an "all-sky" survey of 
rotation measures toward > 10 7 background sources with typical angular 
separations ~ 90", out to redshifts z > 3. 14 Until such information becomes 
available, we cannot rule out the possibility that deflections by intervening 
magnetic fields are significant even for the highest energy CRs. 

3. Candidate astrophysical sources of UHECRs 

A minimum requirement for astrophysical sources of UHECRs is the abil- 
ity to magnetically confine particles of the requisite energies. For particles 
with energy E and charge Z, this implies the condition (i?/pc)(B/lG) > 
(E/10 20 eV)/Z between the system's size R and magnetic field B. Only a 
select few types of objects are known to meet this criterion, among them the 
jets of radio- loud active galactic nuclei (AGNs), gamma-ray bursts (GRBs), 
and clusters of galaxies. 15 Notwithstanding other candidates, in what fol- 
lows, we focus on these three as representative types of potential UHECR 
sources (see Ref. 15 for other possibilities). The actual maximum energy 
attainable under different circumstances must be evaluated by comparing 
the timescales for particle acceleration, usually that for the first order Fermi 
mechanism in shocks, against the timescales for limiting processes such as 
source lifetime, particle escape, adiabatic or radiative energy loss, etc. 

Equally important is the available energy budget. Fig.l shows estimates 
of the kinetic energy output averaged over the universe as a function of red- 
shift z due to AGN jets, GRB explosions and accretion onto clusters, which 
should be proportional to their cosmic ray output. The plotted quantity 
is differential per unit z, dE^/dz = (dt/dz) J L(dn/dL)dL, where L is 
the kinetic luminosity per object and dn/dL is the z-dependent luminosity 
function, with cosmological parameters h=0.7, f2 m =0.3 and flA=0-7. For 
AGN jets, we have made use of the observed radio luminosity function along 
with the observed correlation between the radio and jet kinetic luminosities 
of radio galaxies. 16 GRBs were assumed to occur each with kinetic energy 
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Egkb = 10 53 erg at a rate that follows the star formation history and 
matches the log iV-log S* distribution observed by BATSE 17 (note that this 
estimate is roughly independent of the beaming factor). The three curves 
each for AGNs and GRBs correspond to different evolutionary assumptions 
at the highest z, with only small differences at low z. Cluster accretion will 
be discussed in Sec. 4. 




Fig. 1. Energy budget of candidate UHECR sources: AGN jets (dotted), GRBs (dashed) 
and cluster accretion shocks (solid), the latter separately for each logM as labelled. 



The results at low z can be compared with the observed energy density 
of UHECRs, ~ 3xl(T 20 erg cm" 3 ~ 10 54 erg Mpc" 3 above 10 19 eV. Consid- 
ering further factors for energy loss during propagation, it is apparent that 
whereas AGN jets and cluster accretion shocks have reasonable margins 
to accommodate the energetics of UHECRs, GRBs, with a substantially 
smaller energy budget, require a very high efficiency of energy conversion 
into UHECRs, a fact that has already been noted. 2 ' 18 

For both AGN jets and GRBs, different locations along the outflow can 
be UHECR production sites, since one expects B oc R^ 1 under the naive 
assumption that the ratio of magnetic to kinetic energy is constant. One 
candidate in AGNs is the inner jet region with R ~ 10 16 -10 17 cm and 
B ~ 0.1-1 G, known through observations of blazars to be a site of particle 
acceleration, perhaps due to internal shocks. Estimates show that the max- 
imum proton energy should be limited by photopion interactions with low 
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frequency internal radiation to somewhat below 10 20 eV. 19 Furthermore, 
efheient conversion to neutrons may be necessary to allow the particles to 
escape the jet without suffering adiabatic expansion losses and contribute 
to UHECRs by decaying back to protons outside. A more promising site 
may be the hot spots of powerful radio galaxies, termination shocks where 
large-scale jets are decelerated by the external medium, with R ~ 10 21 cm 
and B ~ 1 mG. Here the maximum energy may exceed 10 20 eV, limited 
by escape. 20 However, note that these particles must further traverse the 
extensive cocoon of shocked, magnetized jet material in order to completely 
escape the system and constitute UHECRs, an issue that has not been ex- 
amined in detail. A further possibility is acceleration by the bow shocks 
being driven into the ambient gas by the expansion of the cocoon. 21 

In order to verify an AGN origin, detailed analysis of observed UHECR 
arrival directions and cross correlations with known source positions will un- 
doubtedly be essential. Nonetheless, given the uncertainties in intervening 
magnetic fields (Sec. 2), it will also be highly desirable to have some means 
to pinpoint individual sources through characteristic, UHECR-induced sig- 
natures of secondary neutral radiation. Neutrinos are ideal in providing an 
unambiguous earmark of high energy hadrons. 22 However, even with km 3 
detector facilities such as IceCube or KM3NeT, it may not be easy to resolve 
individual AGNs. 23 Thus, distinctive electromagnetic signals will also be 
extremely valuable. For radio galaxy hot spots, synchrotron emission from 
UHE protons can produce nonthermal X-rays that could be distinguished 
from other processes such as electron inverse Compton through multiwave- 
length observations. 24 If the medium surrounding the source is sufficiently 
magnetized, UHE protons propagating in the source's vicinity can lead to 
diffuse gamma-ray emission from photomeson-triggered cascades that may 
be detectable by current and upcoming instruments. 25 

Potential locales for UHECR acceleration in GRBs include internal 
shocks, external reverse shocks, and external forward shocks, believed to 
be the emission sites of the prompt X-rays and gamma-rays, optical flash 
and radio flare, and the radio to X-ray afterglow, respectively. 26 The exter- 
nal forward shock could be disfavored due to its ultrarelativistic velocity, 
but this is controversial. 27 For the mildly relativistic internal and exter- 
nal reverse shocks, a different problem is that for the particles to escape 
the acceleration site without significant losses, neutron conversion may be 
required, as with AGN inner jet regions. 28 This entails some reduction in 
efficiency, which cannot be too severe in view of the tight energy budget, 
as discussed above. 
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Identification of high energy neutrino signals from GRBs are facilitated 
through time coincidence, but studying individual bursts in detail may 
be difficult. 29 Thus, photon signatures of UHECR production will also be 
important. 30 We have conducted a detailed study of this issue for internal 
shocks utilizing a comprehensive Monte Carlo code that includes a wide 
variety of relevant processes, and found that in some cases, unique features 
such as synchrotron radiation from muons and protons may be observable 
by GLAST and Cerenkov telescopes (Asano & Inoue, in preparation). 

4. Heavy nuclei from cluster accretion shocks as UHECRs 

In the currently favored picture of hierarchical structure formation in the 
CDM cosmology, all massive clusters of galaxies should be surrounded by 
strong accretion shocks, as a consequence of continuing infall of dark mat- 
ter and baryonic gas. 32 Such shocks should be interesting sites of particle 
acceleration, and have been proposed as sources of UHECRs. 21 ' 33 Here we 
briefly summarize our recent work on this subject invoking UHECR nuclei; 
more details can be found in Rcf. 31. 

For clusters of mass M, the rate of gas kinetic energy dissipation through 
accretion shocks can be estimated as L acc ~ 9 x 1O 45 (M/1O 15 M ) 5 / 3 
erg/s. 34 This can be combined with the Press-Schechter mass function to 
evaluate dE^/dz for clusters of different M as shown in Fig.l. Note that 
due to the hierarchical nature of structure formation together with the non- 
linear nature of gravity, the maximum is reached at z = 0, with ample room 
to supply the UHECR energy budget. 

However, estimates of the maximum energy E maK for protons seem to 
fall short of 10 20 eV by 1-2 orders of magnitude. 33 ' 35 A fiducial cluster 
of M = 2 x 10 15 Mq has shock radius R s ~ 3.2 Mpc and shock velocity 
V s = (4/3)(GM/R s ) 1 / 2 ~ 2200 km/s. The shock magnetic field is taken to 
be B s = 1/iG, as suggested by some recent observations. 36 The timescale for 
shock acceleration is i acc = 20k(E)/V 2 = (20/3)(Ec/ZeB s V 2 ), assuming 
the Bohm limit for the diffusion coefficient k(E) as inferred for supernova 
remnant shocks. 2 To be compared are the energy loss timescales for pho- 
topair and photopion interactions with the CMB, the escape time from the 
acceleration region t csc ~ R 2 /5k(E), and the Hubble time tn- As is clear 
in Fig. 2, for protons -E ma x ~ 10 18 -10 19 eV, confirming previous findings. 

On the other hand, heavy nuclei with higher Z have correspondingly 
shorter t acc , and Fe may be accelerated up to 10 20 eV in the same condi- 
tions, notwithstanding energy losses by photodisintegration with the FIRB 
and CMB (Fig. 2). In order to explore whether nuclei from cluster accretion 
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Fig. 2. Comparison of timescales at cluster accretion shocks for shock acceleration t aC c 
(diagnonal lines), and energy losses from interactions with background radiation fields 
(curves), for protons (thick dotted), He (thin dotted), O (thin solid) and Fe nuclei (thick 
solid). The photopair timescales are denoted separately for p and Fe (dot-dashed). Also 
indicated are the Hubble time tjj (dashed line) and the escape-limited -E ma x (circles). 

shocks can provide a viable picture of UHECR origin, detailed propagation 
calculations of UHE nuclei above 10 19 eV are undertaken, following energy 
losses in the CMB/FIRB and deflections in extragalactic magnetic fields 
(EGMF) for all particles including secondary nuclei arising from photodis- 
integration. We consider EGMF models that trace large-scale structure as 
in Ref. 37, as well as the case of negligible EGMF, although Galactic fields 
are not included. The source density is n s = 2 x 10~ 6 Mpc~ 3 , appropri- 
ate for massive clusters. A fraction Jcr of the accretion luminosity L&cc 
is converted to cosmic rays with energy distributions oc E~ a , and we set 
E max /Z = 10 19 eV, a fair approximation to estimates for each species ob- 
tained by comparing timescales as in Fig. 2. For the elemental composition 
at injection, the He/p ratio is taken to be 0.042. All heavier elements are 
assumed to have the same relative abundances at fixed energy/nucleon as 
that of Galactic cosmic ray sources at GeV energies, 38 and scaled with re- 
spect to protons by the metallicity £ of the accreting gas. We take ( = 0.2 
as suggested by both observations and theory. 39 An additional factor A@ 
for the injected abundance of nuclei with mass number A is introduced to 
take account of possible enhancement of heavier nuclei due to nonlinear 
modification of shock structure by CRs. 40 
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Fig. 3. Observed UHECR spectrum (top) and mean mass composition (bottom) versus 
energy E (1 EeV = 10 18 eV) from cluster accretion shocks for a = 1.7 and fi = 0.5, 
compared with the current data for AGASA 41 (triangles), HiRes 42 (bars) and Auger 43 
(stars). The histograms are the average result over different model realizations for the 
cases with (thick) and without (thin) extragalactic magnetic fields, and the thin curves 
outline the cosmic variance for the former case only. The straight line in the top panel 
denotes a. See Ref. 31 for more details. 

Fig. 3 shows our results for the observed spectrum and composition for 
q = 1.7 and (3 = 0.5, which are quite consistent with the current data for 
HiRes and Auger (and possibly AGASA as well 31 ). Values of a < 2 are nat- 
urally expected at the high energy end for nonlinear shock acceleration that 
accounts for the dynamical back reaction from CRs. 44 The spectral steep- 
ening at > 10 20 eV is due both to propagation losses and the E meix limit 
at the source. Normalization to the observed flux and comparison with the 
available accretion power for M > 1O 15 M0 fixes /cr, which is ~ 0.01 — 0.6 
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for cases with EGMF and ~ 0.004 for the case without. Low values of fen 
may reflect inefficient escape of CRs from the system, which is conceivable 
in view of the converging nature of the accretion flow. CR escape may be 
mediated mainly during episodic merging events that partially disrupt the 
shock and drive outflows of some of the downstream gas. 

The mass composition at < 3 x 10 19 eV is predominantly light and 
consistent with HiRcs claims, 7 while the rapid increase of the average mass 
at higher energies is a clear prediction of the scenario to be tested by the new 
generation experiments. Despite the relative rarity of massive clusters in the 
local universe, strong deflections of the highly charged nuclei in EGMF allow 
consistency with the currently observed global isotropy. On the other hand, 
with a sufficient number of accumulated events, clear anistropies toward a 
small number of individual sources should appear, although this prediction 
is subject to uncertainties in the EGMF and Galactic fields. An aspect of 
this scenario that warrants further study is the spectral domain < 10 19 eV 
and the implications for the Galactic-extragalactic transition region. 18,38 

In this picture, some neutrino production may occur due to confined 
UHE protons interacting with the enhanced FIRB toward the cluster cen- 
ter, 46 along with a component from decay of photodisintegrated neutrons. 47 
Furthermore, we may look forward to very unique signatures in X-rays and 
gamma-rays. Protons accelerated to 10 18 -10 19 eV in cluster accretion shocks 
should efficiently channel energy into pairs of energy 10 15 -10 16 eV through 
interactions with the CMB, which then emit synchrotron radiation peaking 
in hard X-rays and inverse Compton radiation in TeV gamma-rays. Fig. 4 
displays the predicted spectra for a Coma-like cluster, conservatively as- 
suming that UHE proton injection continued only for a dynamical time 
~ 2 Gyr (see Ref. 35 for more details). The detection prospects are very 
promising for Cerenkov telescopes such as HESS, and hard X-ray observa- 
tories such as Suzaku and the future NeXT mission. Photopair production 
by nuclei may also be efficient and induce further interesting signals that 
are worth investigating (Fig. 2). 

Who is the real culprit flinging those nasty UHECRs at us? AGNs, 
GRBs, clusters or some other guy? Through the concerted effort of upcom- 
ing cosmic ray, neutrino, X-ray, gamma-ray as well as radio observations, 
we may finally find out who dunnit! 
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Fig. 4. Spectra of UHE proton-induced photopair emission from the accretion shock of 
a Coma-like cluster, for B s =0.1, 0.3 and 1 jiG. The sensitivities for a 1 degree extended 
source are overlayed for HESS, GLAST, Suzaku XIS+HXD, and NeXT HXI+SGD. 
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